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abstract 


la  this  paper  we  analyze  Che  numerical  aspects  of  the  various  methods 
Chat  have  been  used  to  analyze  chin  wire  antenna  problems-.  First  we  derive 
some  properties  of  the  thin  wire  E-field  Integral  operator.  Based  on  those 
properties  we  unify  the  various  iterative  methods  used  to  find  current  distri¬ 
bution  on  thin  wire  structures.  An  attempt  has  been  made  to  resolve  the  ques¬ 
tion  of  numerical  stability  associated  with  various  entire  domain  and  subdo¬ 
main  expansion  functions  in  Galerkin's  method.  It  has  been  shown  chat  the  se¬ 
quence  of  solutions  generated  by  the  Iterative  methods  monotonlcally  approaches 
Che  exact  solution  provided  the  excitation  chosen  for  these  problems  are  in 
the  range  of  E-field  operator.  Such  a  statement  does  not  hold  for  Galerkln's 
method  since  the  inverse  operator  is  unbounded.  Moreover  if  the  excitation 
function  is  not  in  the  range  of  the  operator  the  sequence  of  solutions  form  an 
asymptotic  series.  Examples  have  been  presented  to  illustrate  this  point. 
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I.  INTRODUCTION 


Over  the  past  thirty  years  several  methods  have  been  developed  by 
various  researchers  to  analyze  scattering  and  radiation  from  thin  wire 
structures.  In  this  presentation  we  investigate  the  properties  of  the 
integrodifterential  equations  by  various  techniques  which  impose  boundary 
conditions  on  the  electric  field  only.  Se  we  restrict  our  discussions  to 
the  E-fleld  integral  equation. 

For  simplicity  let  us  focus  our  attention  on  the  reradiation  of 
electromagnetic  waves  by  a  thin  wire  of  radius  a  and  of  length  L  centered 
at  z  »  0.  In  this  analysis  we  will  assume  that  the  antenna  is  quite  thin 
and  hence  the  circumferential  variation  of  the  current  on  the  antenna 
can  be  neglected.  He  shall  also  neglect  any  $  component  of  the  current 
on  the  structure.  Furthermore,  there  are  no  internal  resonances  of  the 
thin  wire  at  the  frequency  of  operation.  With  these  assumptions  we 
derive  the  properties  of  the  E-fleld  operator  for  both  the  Pocklington 
and  Hallen’s  equations. 


2.  PROPERTY  OF  THE  POCKLINGTON  E-FIELD  OPERATOR 


By  assuming  a  time  variation  of  the  form  exp  (jut)  the  Pocklington 
Integral  equation  for  the  current  on  the  surface  of  antenna  can  be  written 
as  (1) 

+L/2 


2ir 


+L/2  2n 


k2  /  dz  /  d<f>  I(z')  G(z,z’)  +  — 2  /  dz'  fd<p  I(z')  G(z,z')  “  -jutfireE^z) 


-L/2  o 


dz  -L/2  o 


for  -  j  <  z  <  +  y 


where  G(z,z’)  -  ^ 


(1) 

(2) 


R  -  /“!  ,.2i(  2  .24 

(z-z')  +  4a  sin 


(3) 


E  (z)  »  incident  field  on  the  antenna  (4) 

In  the  operator  form  (1)  can  be  written  as 

PI  A  k2  Px  I  +  P2  I  -  V  (5) 

We  next  investigate  if  there  exists  a  constant  C  Independent  of  I  such 
that  [2,  p.296] 
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||PI||  <C 


I  -  1 


(6) 


If  such  a  constant  C  exists  which  is  the  0f  all  possible  j  [plj  | 

with  the  constraint  | |l| |  >1  then  we  say  the  operator  is  bounded  with 

respect  to  that  norm  |  ( *  1 1 .  The  norms  that  we  shall  be  dealing  with  are 

2  2 
the  if  norm  and  the  Chebyshev  norm  .  The  norm  is  defined  aa 


l|l|L-  I  V  |K*>|2 

it1  -L/2 


(7) 
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and  the  Chebyshev  norm  is  defined  as 

IUII  -  L  |I(*)|  (8) 

T  L  <  z  <  +  j 
~  2 

When  no  subscripts  are  used  either  nora  is  implied.  If 

we  are  using  the  &  norm  then  we  are  restricting  the  domain  of  the  oper- 

2 

ator  P  to  elements  which  are  in  (or  square  integrable) .  This  does  not 
imply  that  I(z')  cannot  be  infinite  within  the  range  -  y  r '  £  +  j  • 

However  only  those  type  of  singularities  are  permitted  in  I(z')  which  are 
square  integrable.  Any  function  which  is  not  square  integrable  is  ex¬ 
cluded  from  the  domain  of  P  (as  they  are  not  in  «^).  On  the  other  hand, 
if  we  use  the  Chebyshev  norm  then  the  function  has  to  be  bounded.  Under 
the  Chebyshev  norm  any  unbounded  function  cannot  be  in  the  domain  of  the 
operator.  Thus  the  function  log  z  is  in^  but  not  in  the  domain  of  func¬ 
tions  bounded  under  the  Chebyshev  norm.  Physically  then  convergence  of  a 
sequence  of  functions  under  the  norm  yields  least  squares  convergence 

whereas  convergence  under  the  Chebyshev  norm  yields  pointwlse  convergence. 

Examination  of  the  first  part  of  the  Integral  in  (1)  reveals  that 
the  kernel  has  a  singularity  which  can  be  observed  by  rewriting  the  kernel 
as  [3,  p.  141] 

2ir  2ir 

/  d*  G(z,z')  -  J  <  ~JkR?  d*  -  log  |z  -  *•  |  +  G.  (9) 

o  o  * 

where  contains  terms  which  are  bounded  and  hence  square  Integrable. 

The  singularity  of  the  kernel  in  P^  is  manifested  through  the  function 
log  |s  -  «' |.  Since  a  lograthmlc  function  is  square  integrable  we  find 
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+L/2  +L/2 

-In  (i««  l-’l  ♦o,)V/2-c< 
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where  C  is  a  constant.  Hence  is  bounded  under  the  norm.  Furthermore 
P^  is  a  Hilbert  Schmidt  operator  as  it  has  a  square  integrable  kernel 
[2,  p.  352].  It  can  also  be  shorn  that  a  Hilbert  Schmidt  operator  is 
a  compact  operator  [2,  p.  353].  Under  the  Chebyshev  norm 


L  .  L 

T  -  Z  -  +  2 
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+L/2 

/  dz'  log  |2-z’| 

1  +"“ 

1  L  T 
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-L/2 

i  2  —  —  2 
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|  (z+  j)  log(z  +  y)  +  (  j  -  a)  log  (  Y  -  z)  -  L|  +  constant 


_<  M  (a  constant)  for  all  -  j  <  z  <  +  |  (11) 

Hence  the  operator  P^  is  also  bounded  under  the  Chebyshev  norm, 
lest  consider  the  second  integral  in  (1).  He  have 
+L/2  -H./2 

1  •  h  [£'  &  (1««  I—' I  +  °2‘  •  &  h  +  o2> 

_ (12) 

where  the  bar  over  the  aecond  Integral  represents  a  principal  value.  It 

is  clear  that  the  operator  Pj  in  (12)  is  unbounded  under  the  Chebyshev 
31  L 

norm,  because  is  unbounded  as  z  ■*  ♦  r  ,  Thus  there  exist  no 

constant  C  for  all  -  j  <  s  <  ♦  j  such  that 

IIp2Nt lc  <m 


<■»  2  3X 

Also  the  operator  Is  unbounded  under  the  ^  norm  as  ,  and 

J 

-fa  {log  |z  -  zf|  +  G^}  are  not  square  lntegrable. 

31 

However  If  the  antenna  has  no  edges  (i.e.  g^is  everywhere  bounded 

and  square  lntegrable)  then  it  can  be  shown  that  the  operator  is  bounded 
2 

under  the  X  norm.  Sneddon  has  shown  through  Theorem  8  [10,  p.  234]  that 
if  -f(z)  is  square  integrable  over  -  ^  £  z  y  and  zero  everywhere  else 

then  the  formula  fH(z)  *  -  1  j  f  (t)  log  j~|Zl  dt 

defines  almost  everywhere  a  function  J„(z)  which  is  also  square  integrable,  and 

n 


I l/(z)l 


IV2) 


Hence  if  the  antenna  has  ’lo  sharp  edges  then  the  operator  is  bounded 
2 

under  theX  norm  as  illustrated  by  the  above  theorem.  In  other  words. 


if  the  antenna  has  end  caps  then  the  Pocklington  E-field  operator  is  bounded. 


Hence  the  Pocklington  E-field  operator,  P  »  P^  +  P^  is  unbounded  if 
the  scatterer  has  a  sharp  edge  at  which  the  charge  density  becomes 
infinite.  However  if  the  antenna  structure  is  quite  smooth,  so  that  at 


each  point  the  tangent  to  the  surface  can  be  defined  uniquely,  then  the 
Pocklington  E-field  operator  P  -  P^  +  P^  is  bounded  under  the  norm. 
Next  we  consider  the  Hallen  E-field  operator. 


3.  PROPERTY  OF  THE  HALLEN  E-FIELD  OPERATOR 

Hallen  transformed  Pocklington' s  equation  as  given  by  (1)  into  the 
following  integral  equation 
+L/2  27V 

/  dz'  /  d$  I(z')G(z,z')  •  D  cos  kz  +  F  sin  kz  + 

-L/2  0 

z 

^  f  E*(z)  sin  k(z  -  z')  dz'  (14) 

“  -L/2 


-3- 


where  D  and  F  are  constants  which  depend  on  how  the  current  goes  to  zero 
at  the  ends  of  the  antenna.  We  define  the  Kallen  E-field  operator  as 


+L/2  2ir 

HI  -  /  dz'  f  dip  I(z')  G(z,  z')  -  P.  I  [from  (1)1  (15) 

-L/2  o  1 

H  ence  the  operator  H  is  bounded  both  under  the  and  the  Chebyshev  norm. 
Also  H  is  a  compact  operator  under  thej£  norm. 

It  is  Important  to  note  however,  that  the  unknown  I(z')  is  hidden  in 
the  constants  D  and  F.  To  illustrate  this  further,  if  we  consider  a  delta 
gap  excitation  for  the  antenna  then  (14)  becomes  [4,  p.  325] 

+L/2  2ir 

/  dz'  /  d<J>  I(z’)  G(z,  z’)  ■  A  sin  k  |z|  +  D  cos  kz 
-L/2  o 

where  A  is  a  known  constant  and  0  is  unkown.  Observe  at  z  ■  0 
+L/2  2ir 

D  /  dz’  /  dp  I(z')  r.(o,  z') 

-L/2  O 

If  the  operator  H  is  bounded  so  will  be  the  constant  D.  If  one  wishes 
then  perhaps  one  can  transfer  D  to  the  left  hand  side  of  the  equation 
and  thus  form  an  additional  part  of  the  operator  H.  But  sine  e  D  is  a 
part  of  the  operator  H,  whatever  bounds  hold  for  H  also  hold  for  0. 

Next  we  estimate  a  bound  for  | |h| |  both  under  the  5^  and  Chebyshev 

norm. 


He  observe 
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■  max  [2  log  log  ~  ]  +  @(a^)  [3,  p.  144] 

fl  3 

2  2  L 

where  9  (a  )  denotes  terms  of  the  order  of  a  .  Now  if  —  »  1  then 

.  a 

INI  <2  log  t 
T 


Under  the  ^5  norm 


+L/2  +L/2  2ir 

lHl  L2  1  I  /  dz  /  dz'  t  /  d<f 

•C*  -L/2  -L/2  o 


exEizj 

2ttR 


.  2  1/2 
L>  ] 


+L/2  +L/2  2ir  2  1/2 

<  [  /  dz  /  dz'  {  /  d<J>  }  1 

-L/2  -L/2  o  Z1TR 


<  ®  [  —  log  — ]  +  other  terms 
a  a 


The  above  Integrals  may  be  evaluated  numerically  to  give  a  more  accurate 
estimate  for  IN^g* 

As  H  is  bounded  operator  unlike  P,  it  may  be  computationally  much 
easier  to  solve  Hallen  integral  equation  than  Pocklington  integral 
equation. 


4.  SOLUTION  OF  HALLEN  E-FIELD  INTEGRAL  EQUATION 
4.1.  By  Iterative  Methods 

It  is  well  known  that  if  H  is  a  compact  invertible  operator  (under 

2 

norm)  on  an  infinite  dimensional  space  then  its  inverse  is  often 
unbounded  [3,  p.  353].  Hence  the  problem  of  the  solution  of  (14)  in  the 


-7- 


if  norm  is  ill-posed*  If  a  problem  is  ill-posed  under  the  norm  then 
it  is  definitely  ill-posed  under  the  Chebyshev  norm.  However  it  can 
be  regularized  in  the  following  way.  We  take  (14)  and  cast  it  in  the 
form 

HI  «  Q  (19) 

and  generate  the  sequence  {i^} 

I  -  [U  -  x  H]  I  .  +  T  Q  (20) 

n  n— l 

with  a  starting  guess  of  Iq  3  Q  and  U  is  the  identity  operator.  The 

sequence  I  generated  by  (20)  converges  to  a  solution  I  which  satisfies 
n 

HI  ■  Q  for  all  Q  in  the  range  of  H  [1,  p.  196].  The  sequence  generated 
by  (20)  always  converge  to  I  provided 


-8- 


(when  convergence  is  desired  in  ^  norm) ;  and 

,  M  h|  L 


or 


2  log  — 
a 


(23) 


i  w  i  *  rif 

(when  convergence  is  desired  in  Chebyshev  norm) . 

Observe  if  solution  is  desired  and  one  chosses  t  according  to  (23) 

it  is  quite  possible  that  the  iterative  method  described  by  (20)  may 

not  converge  at  all.  Therefore  if  we  are  solving  for  the  current  distribution 

on  a  thin  wire  due  to  a  delta  excitation  the  value  of  x  chosen  should  be 

0  [—  log  — ]  instead  of  the  popularly  used  value  of  T  =  2  log  — .  This  is 
a  a  a 

because  the  current  has  a  logarithmic  singularity  at  z  =  0  and  so  convergence 

2 

of  the  various  iterations  can  be  guaranteed  only  in  an  sense. 

Hallen  in  his  classic  iterative  scheme  chose  the  value  of  T  as  given 
by  (23)  [3,  4],  A  detailed  description  on  how  D  is  solved  for  at  each 
iteration  is  described  in  detail  [4,  p.  326],  Observe  that  if  ^  ^ 

or  | |U  -  x  H| |  <1  then  the  iterations  defined  by  {(8,149)  and  (8.150)  of 
[4]}  would  always  converge  for  any  starting  guess  1^. 

Other  researchers  have  chosen  different  values  of  x.  For  example 
Gray  [5]  chose 


-  -  Real  [2  log  -  -  2y  -  2  log  —  -  jir  +  2E1  (-  ^ 
x  a  l 


)] 


where  y  is  Euler’s  constant  and  Ei  is  the  exponential  integral, 

and  Middleton  [6]  decided  to  make 
L/2 

^  ■  /  G(L  -  ~  f  z’)  sin  kzf  dzf. 

0 

Whereas  Siegel  and  Labus  [7]  chose 

—  -  2  log  -  -  C  in  (kL)  -  1  - 
T  a  RL 


(24) 

King 

(24. a) 


(25) 
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Where  Cin  is  the  special  form  of  the  cosine  integral.  Finally  Schelkunoff 
[3]  after  a  careful  analysis  decided 

-  -  2  log  \  -  Cin(kL)  -  1  -  -^kL  -j  Sin(kL)  +  j  kL(25.a) 

In  general,  it  really  does  not  make  any  difference  whatsoever,  what 

value  of  T  one  chooses,  one  is  guaranteed  to  have  pointwise  convergence 

or  convergence  in  the  mean  depending  on  whether  one  chooses  T  according 

to  (23) or  (22).  This  of  course  assumes  that  a  solution  to  the  problem 

exist,  i.e.  Q  is  in  the  range  of  H. 

In  summary,  the  iterative  method  convert  HI  =  Q,  a  Fredholm  equation 

of  the  first  kind  to  I  =  B  I  _  +  tQ,  as  Fredholm  equation  of  the 

n  n—  I 

second  kind.  The  advantage  of  the  equation  of  the  second  kind  is  that 
1 1  [U  -  B]  |  |  is  not  only  bounded  but  its  inverse  \  | [\3  -  B]  |is  also 
bounded,  provided  unity  is  not  an  eigenvalue  of  B.  Mathematically  one 
has  regularized  the  problem  by  the  introduction  of  the  parameter  t. 

With  this  regularization  scheme  the  convergence  of  the  sequence  1^  is 
mono tonic. 

Finally,  we  conclude  by  noting  that  as  the  iterative  process 
continues  the  unknown  constants  D  and  F  in  (14)  are  determined  as  out¬ 
lined  in  [4]. 

4.2#  By  Galerkin's  Method 

The  next  generation  of  the  methods  were  developed  primarily  by 
Harrington  [8]  under  the  generic  name  of  "moment  methods."  This  very 
popular  versatile  method  has  been  excellently  documented  in  [8].  In 
Galerkin's  method,  the  unknown  function  I  is  expressed  as 
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N 

Ij,(z)  “  l  ^(z)  (2d 

where  f^(z)  are  known  functions  which  may  extend  from  -  j  z<_  +  |  or 

could  span  only  a  partial  portion  of  the  domain  of  z,  i.e. 

-  \  <  ai  —  Z—  °2  <  +  \  #  ***  ^ormer  case  ^  ’s  become  entire  domain 

functions  whereas  in  the  latter  Y  *s  are  called  sub  domain  basis  functions. 

We  solve  for  I  (z)  by  solving  for  the  unknowns  a.  in  (26) .  We  also 
n  i 

convert  the  Infinite  dimensional  problem  HI  ■  Q  to  a  finite  dimensional 

problem  by  replacing  I  with  1^  ,  i.e.  we  solve  the  following  equation 
N 

£  a  HY  ■  Q  in  the  finite  dimensional  space  spanned  by  the  basis 
i-1  1  1 

functions  V,,  i  ■  1,  2,  ...t  K.  We  next  find  a  unique  solution  in  finite 

N 

dimensional  space  by  weighting  the  residual  J  a.  -  Q  to  zero  in 

i*l  1  1 

the  following  way 
N 

I  a .  <1^.,  f,>  -  <Q,  ¥,  >  for  j  -  1,  2 . N  (27) 

i«l  1  3  J 


In  a  matrix  form 


[G]  [a]  -  [V]  (28) 

where  [G]  -  [<1^,  4^>] 

[V]  -  [<Q, 

and  the  inner  product  is  defined  as 
+L/2 

■  /  dz  4».  (z)  ♦.  (z) 

1  J  -L/2  1  3 

The  unknown  a  's  in  (28)  are  obtained  as 

[a)  -  IG]"1  [Z]  (29) 
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The  next  question  that  normally  arises  in  whether  the  sequence  defined 
in  (26)  approaches  any  limit  I  as  N  -►  °°.  And  secondly  whether  I  satisfies 
the  equation  HI  ■  Q.  We  cannot  talk  about  convergence  in  the  Chebyshev 
metric  [as  defined  in  (8)]  because  a  Chebyshev  norm  cannot  be  derived 
from  an  inner  product  [2.  p.  272].  In  other  words  in  an  inner  product 
space  we  cannot  define  a  Chebyshev  norm.  Hence  we  shall  be  talking  about 
only  the  norm  for  Galerkin's  method.  So  we  shall  be  discussing  about 
convergence  in  the  mean.  Galerkia's  method  guarantees  the  weak  conver¬ 
gence  of  the  residuals  [from  (27) ]9  i.e. 


Lt 

N  -►  « 


-  Q>  y  * 0 


for  j  *  1,  2,  .... ,  N 


(30) 


However  if  H  is  a  bounded  operator  (i.e.  | |h| |  7 

*£r 

(30)  implies  strong  convergence  of  the  residuals. 


<  a  constant  <  <*>)  then 
i.e. 


1 1“»  -  OlJj.1  * 0  <3l) 

This  has  been  proved  by  Mikhlin  [9].  Physically,  (31)  implies  that  as 
N  the  tangential  electric  field  on  the  surface  of  the  conductor 

converges  to  zero  in  a  least  squares  fashion. 

Unfortunately  in  Galerkin's  method  the  convergence  of  the  residuals 
to  zero  in  (31)  does  not  imply  the  convergence  of  1^  to  a  solution  1  of 
HI  “  Q.  The  convergence  of  1^  +  I  in  the  domain  of  H  is  possible  if  and 
only  if  | |H_1|  2  is  bounded,  as 

IIS  -  i  U®"1!^  •  U”n-<>li.2  <32> 

So  it  1 1 H  3|L2  1.  unbounded,  mo  though  tho  rosldunls  go  to  o.ro ,  tho 
sequence  of  solutions  1^  may  not  converge  to  1.  This  is  in  contrast  to 
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'the  iterative  methods  where  mono  tonic  convergence  to  I  is  guaranteed  if 
T  and  Q  are  chosen  as  prescribed. 

Since  | |H  |  -  *8  unbounded  in  this  case,  the  application  of 
qCz 

Galer kin’s  method  to  HI  *  Q  may  not  guarantee  that  |  |l^  -  l|  4  0  as  N  + 
In  other  words,  there  is  no  quantitative  way  to  describe  the  convergence 
of  1^  -►  I  as  various  expansion  functions  are  chosen  for  4^.  Hence  we 
address  the  question:  For  a  fixed  order  of  approximation  N,  h  w  should 
one  choose  a  set  of  expansion  functions  V  such  that  the  round-off  and 
the  truncation  error  in  the  numerical  computation  of  a  in  (29)  is  a 
minimum? 

Suppose  the  Gram  matrix  £  is  generated  by  the  basis  functions 
[Ey  *  <4^,  4*  >]  then  we  show  in  the  appendix  (Theorem  3)  that 

/V 

cond  [G]  £  cond  [H],  cond  [E]  (33) 

l.e.  the  condition  number  of  the  Galerkin  matrix  G  in  (28)  is  bounded 
by  the  condition  number  of  the  operator  H  in  the  finite  N  dimensional 
space  and  the  Gram  matrix  E.  Equation  (33)  is  valid  only  in  the  finite 
N  dimensional  space  spanned  by  ¥  .  It  is  important  to  note  that  even 
though  H  may  not  have  any  eigenvalues  in  an  infinite  dimensional  space, 
it  has  at  least  an  eigenvalue  on  a  finite  dimensional  space  [2,  p.  332]. 

If  the  homogeneous  equation  HI  ■  0  has  only  the  trivial  solution  I  ■  0  and 
| | H | |  2  is  bounded  then  cond  [H]  <  °°  and  the  inequality  in  (33)  has  some 

meaning  because  the  right  hand  side  of  (33)  can  never  be  infinity. 

So  (33)  directly  implies  the  following:  i)  Use  of  an  orthonormal 


set  of  basis  functions  for  the  current  Implies 

A 

cond  [G]  £  cond  [H] 


(34) 
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i.e.  the  problem  would  not  be  worse  conditioned  as  the  original  problem. 

For  this  case  cond  [E]  ■  1.  This  will  happen  for  subdomain  basis  functions 

like  pulses  or  entire  domain  orthonormal  basis  functions  like  /~2  Sin  m  z 

ir 

for  m  •  1,  2,  .  N  .  (34)  also  implies  that  the  solution  of  HI  ■  Q  by 

Galerkln's  method  in  a  finite  dimensional  space  may  be  a  better  conditioned 
problem  than  the  original  problem  posed  in  the  finite  dimensional  space  N. 
This  definitely  should  be  a  very  strong  point  for  Galerkln's  method. 

Also  from  (34)  there  is  no  way  to  tell  whether  the  Galerkin  matrix 
G  associated  with  the  entire  domain  basis  functions  would  be  more  ill- 
conditioned  than  the  Galerkin  matrix  associated  with  the  pulse  functions. 

11)  Use  of  subdomain  basis  functions  like  triangles  or  piecewise 
sinusoids  may  deteriorate  the  condition  number  of  the  Galerkin  matrix 
[G]  than  that  of  the  original  problem.  This  is  because  cond  [E]  >1  for 
these  cases. 

For  the  case  when  's  are  chosen  as  piecewise  triangles,  then  E 
is  a  tridiagonal  matrix  of  the  form 

P  Q  0  . 

QPQ. 

0  Q  P  . 

-  . 

where  P  ■  and  Q  ■  and  Az  «  .  Since  the  jth  eigenvalue 

of  a  tridiagonal  matrix  is  given  by  (1.  p.  70] 

•  P  +  2Q  cos  (JJj) 

we  have 

cond  [EJ  <  1PI  ± 2  W  -  3 

triangles"  |P|  -  2  |q| 
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(36) 


Bence  for  all  dlaenslon  N  the  Galerkin  matrix  due  to  piecewise  triangle 
expansion  functions  may  have  a  condition  number  which  at  most  can  be  three 
times  as  that  of  the  original  problem,  l.e. 

cond  [G]  <  3  cond  [H] 


For  the  piecewise  sinusoids  however, 


_  _  2k  As  -  Sin  2k  Az  ,  »  _  Sin  k  A  z  -  k  A  z  cos  k  A  z 
r  *  - s -  and  Q  »  - -  ~  -  -  ■  - 

2k  Sin  k  A  z  2k  Sin  k  A  z 


(37) 


where  k 


2ir 


wavelength 


.  In  this  case  cond  [E]  is  bounded  by 


Cond  [E]  <_ 

Sinusoids 


1 2k  A  z  -  Sin  2k  A  z|  +2  | Sin  k  A  z  -  k  A  z  cos  k  A  z| 
|2k  A  z  -  Sin  2k  A  z|  -  2  |Sln  k  A  z  -  k  A  z  cos  k  A  z| 


(38) 


In  the  limit  A  z  •*  0 


cond  [Ej  <  3  (39) 

Sinusoids 

Thus  (39)  implies  that  as  the  dimension  of  the  problem  becomes  large  the 
Galerkin  matrix  due  to  piecewise  sinusoids  are  no  less  numerically  ill- 
conditioned  than  the  matrix  produced  by  piecewise  trlsngles.  It  may  be 
quite  possible  that  for  a  particular  value  of  K  the  Galerkin  matrix  due 
to  piecewlae  ainusoldal  functions  may  be  better  conditioned  than  that  of 
the  piecewise  triangles  or  even  than  that  of  the  pulse  functions. 

In  the  above  analysis  an  attempt  has  been  made  to  provide  a  worst 
case  theoretical  bound  for  tha  condition  number  of  the  various  matrices 
of  interest. 

It  is  important  to  stress  that  the  problem  we  have  addressed  here  is 
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not  which  set  of  basis  functions  would  provide  the  best  approximation 
for  the  current,  but  which  type  of  expansion  functions  would  give  rise  to 
a  well  conditioned  Galerkin  matrix  G  which  will  be  easy  to  invert  numerically. 
This  is  because  truncation  and  round-off  error  associated  with  the  solution 
of  (28)  is  directly  related  to  cond  [G]. 

5.  IS  A  SOLUTION  POSSIBLE  IF  THE  EXCITATION  IS  NOT  IN  THE  RANGE  OF  THE 
OPERATOR? 

He  discuss  the  question  of  existence  of  a  solution  for  the  current 
on  the  antenna  structure  when  we  try  to  excite  it  with  a  source  which  is 
not  In  the  range  of  the  operator  H.  Clearly,  if  the  excitation  is  not  in 
the  range  of  the  operator  then  mathematically  a  solution  does  not  exist. 

But  numerically  one  could  always  find  a  solution  to  the  integral  equation. 

This  numerical  solution  has  some  very  interesting  properties  as  outlined 
in  Theorems  4  and  5  (in  the  appendix) .  If  we  try  to  numerically  solve 
an  integral  equation  HI  »  Q  with  Range  of  H,  then  the  sequence  of 
solutions  1^  diverges  even  though  the  residuals  HI^  -  Q  associated  with 
HI  “  Q  may  approach  zero  monotonically.  This  has  been  proved  in  Theorem  4 
of  the  appendix.  In  theorem  5,  we  develop  further  properties  of  the 
solution  Ijj.  There  we  prove  that  the  sequence  1^  indeed  form  an  asymptotic 
aeries.  The  asymptotic  series  has  the  property  that  it  converges  at 
first  and  then  as  more  and  more  tens  are  included  in  the  series,  the 
series  actually  diverges.  Even  though  the  theorems  4  and  5  have  been 
proved  for  the  iterative  methods,  they  are  alao  valid  for  Galerkin* s 
method.  He  now  present  some  examples  to  illustrate  when  Q  is  in  the 
range  of  operator  and  when  it  is  not. 

•  flr«t  example  consider  the  radiation  problem  where  an  antenna 
of  length  L  and  radius  a  is  excited  by  a  delta  gap  at  the  center.  The 
corresponding  Ha lien’s  integral  equation  has  the  following  form  [4,  p.  3211 


where  A  is  a  known  constant  and  D  Is  unknown.  It  has  been  shown  by  Wu 
[4,  p.  322]  that  the  solution  I(z)  of  the  above  equation  has  a  logarithmic 
singularity  at  z  »  0  given  by 

I(r)  “  '  J  12ot  V  U°8  ]7[+  ••**1 

Hence  I(z)  to  the  delta  excitation  is  square  integrable  but  is  not  bounded. 

So  the  delta  function  excitation  would  be  in  the  range  of  the  operator 
2 

if  the  AT  norm  is  used.  The  delta  function  excitation  would  not  be  in 
the  range  of  the  operator  if  the  Chebyshev  norm  is  used. 

As  a  second  example  consider  the  scattering  from  a  wire  of  length  L 
and  radius  a  irradiated  by  a  broadside  incident  plane  wave.  In  this 
case  the  integral  equation  would  be  of  the  form 

+L/2  2ir 

/  dz'  /  d$  I(z')  G  (z,z')  “  D  cos  k  z  +  F  Sin  k  z  +  Y 
-L/2  o 

where  Dt  F  and  T  are  constants.  By  superposition  principle  there  must 

exist  a  current  Ly(z' )  such  that 
+L/2  2ir 

/  dz'  /  d$  ly(z')  G(z,z')  “  a  constant 
-L/2  o 

+L/2  2lt  dL(*') 

/  dz'  /  d*  ~ jp-  G(z.z')  -  0  (40) 

-L/2  0 

If  we  assume  that  the  function  I^(z')  is  square  integrable  then  we 
can  expand 

R 

Iyt*')  -  I  \  *k<*)  (41) 
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whare  ?k(z)  -  1  for  *k  <  *  <  zfc+1 
•  0  otherwise 

By  substituting  (41)  into  (40)  yields 
H  +L/2  2it 

l  \  I  dz’  /  d«  (G(z,  zfc)  -  G(z,  zk+1)]  -  0 
k-1  -L/2  0 

2 

Since  G(z,zk)  1*  G^.z^j)  t*len  °k  *  ®  for  ^  "  1»  2,  .....  N.  Hence  an  & 
solution  does  not  exist  for  this  problem.  However  if  we  impose  the 
additional  constraints  I(+  ■  0  then  0  cos  kz  +  F  Sin  kz  +  Y  may  become 

an  element  in  the  range  of  the  operator. 

Whether  in  fact  a  forcing  function  is  in  the  range  of  the  operator 
is  difficult  to  verify  both  theoretically  and  numerically.  If  a  forcing 
function  is  not  in  the  range  of  the  operator,  theoretically  we  should 
obtain  a  solution  which  diverges  in  an  asymptotic  sense.  However  this 
postulate  may  be  difficult  to  verify  numerically  for  certain  problems. 

As  an  example  consider  the  partial  sum  of  the  series 


-  1  + 


+ 


1 

N 


The  partial  sum  S^  diverges  u  S  +  *.  This  is  because  if  we  look  at  the 
following  M  terms  of  the  series  we  find 


if  we  program  the  series  on  the  computer  sad  ask  the  computer  to 


give  us  a  result  when  the  addition  of  the  N+l  term  does  not  change  the 
partial  sum  by  10  (say)  we  would  get  a  convergent  result! 

In  conclusion,  we  must  try  to  learn  theoretically,  as  much  about  the 
problem  as  possible.  Numerical  methods  may  be  applied  as  a  last  resort 
as  it  may  be  the  only  way  to  obtain  a  solution  easily.  The  convergence 
of  the  numerically  computed  results  is  determined  to  a  large  extent  by 
the  theoretical  analysis  of  the  problem  rather  than  apparent  convergences 
in  numerical  computations. 

6.  CONCLUSION 

In  summary,  we  have  brought  out  the  following  features. 

1)  The  thin  wire  E-field  Pocklington  Integral  operator  is  unbounded 
whereas  the  Hallen  E-field  operator  is  bounded. 

2)  The  inverse  operator  for  Hallen  E-field  integral  equation  is  unbounded. 

3)  K  unified  theory  for  the  various  iterative  methods  showing  how 
the  Fredholm  equation  of  the  first  kind  has  been  converted  to  a  Fredholm 
equation  of  the  second  kind  is  presented. 

4)  The  conditions  under  which  the  iterative  methods  converge  both 

2 

for  the  &  norm  and  the  Chebyshev  norm  has  been  presented. 

5)  The  monotonic  rate  of  convergence  of  the  sequence  of  solutions 
associated  with  iterative  methods  have  been  established  for  certain 
values  of  t  and  for  Q  e  R(H). 

6)  The  numerical  stability  in  the  solution  of  the  matrix  equations 
for  Galerklnfs  method  for  various  expansion  functions  is  examined  and 

7)  The  sequence  of  solutions  1^  forms  an  asymptotic  series  for  both 
the  Iterative  and  Galerkln's  method  when  Q  t  R(H) . 
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8.  APPENDIX  v 
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and  so 


from  which  it  follows  that 

A 

X  [G]  <  I |H| I •  X  [E] 
max  —  1 1  *  1  max 


Also  since 


| <HI,I> |  > 


1 1  i 


<EI,  I>  >  Amin[E]*l 

l|H"1||~  I  I H-1 1  I 


1  <HI,  I>|  >  Xmln[E] 

<I,I>  “  I  I H~ 1 1 


from  which  it  follows 


WG>  > 


Hence  we  have 


X  [C] 

max 


Xmin[Gl 


*  a  ,  X  [E] 

<  1 1 H f  |  •  |  |H  X|  |  - 


WEJ 


cond[G]  _<  cond[H] •  cond  [E]. 

Theorem  4 :  If  Q  £  R(H)  then  the  sequence  of  approximations  I  generated 
by 

I  . ,  -  BI  +  Q 
n+1  n 

with  the  initialization  1^  *  Q  yields  the  following  relationships 

i)  lim  | | R  -  rI |  -  0,  where  R  -  HI  -  Q 
Mm  nTX  n  n  n 


ii)  lim  III 


m  CO 


Proof ;  i)  We  have 


R  ..  -  R  -  H[I  .  -  I  ] 

n+1  n  n+1  n 

and  since  the  operator  H  is  bounded  [i.e.  |  j  H | |  <  M]  we  have 


ll*+,  -  *11  <  «  Hi  «  -  i  Il3<  Ik 


n+1  “  Gnl 


<  M-||B  -  U||  •  ||en||  <  M.  | I B  -  U | I •  {|!b||}-||£c 
Hence  lim  | -  Rni I  =  0  as  lim  {| |Bn| |}  -►  0. 


n-*» 


n-*» 


ii)  If  lim  |  1 1  ||  =  00  does  not  hold,  then  we  have  lim  |  |l  |  |  <  » 
nr*»  n  n n 

(i.e  a  bounded  sequence).  Thus  there  is  a  subsequence  I1  which  is 

n 

bounded  in  norm.  Now  if  we  put  the  operator  equation  in  a  Hilbert  space 

2 

setting  (now  we  can  only  talk  about  the  norm  only)  and  since  a 

Hilbert  space  is  weakly  compact  [2]  one  can  always  extract  from  1^  another 

sequence  1^  which  converges  weakly  to  some  element  I  of  the  Hilbert 

space,  i.e.  I”  ^  I.  Also  we  have  from  i)  lim  HI  ^  Q  (strong 

ir*00  n 

convergence  in  norm) .  However  as  H  is  a  bounded  operator  we  have 

lim  H  I”  $  HI 
d-n»  n 

and  also  ^ 

11a  H  I  +  Q 
n 

Since  the  weak  and  strong  limits  of  a  sequence  must  concide,  H  I"  *  Q. 

n 

This  means  Q  £  R(H),  a  contradiction. 

Theorem  5:  The  sequence  I  as  derived  in  ii)  of  Theorem  4  indeed 
forms  an  asymptotic  series  (i.e.  for  good  accuracy  the  series  has  to  be 
truncated  after  a  finite  number  of  terms  otherwise  the  results  may  be 
worse) . 
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Proof:  To  demonstrate  the  source  of  divergence  in  I  we  assume 
- -  n 


HI  -  Q  with  I  -  Q  +  AQ 
o  o  1  o 


I  -  I  -  [B]  I  +  0  .  for  n  <  1 

n  o  n-1  — 


where 


0 _  ,  -  1  [B]1  AQ 


If  Q  +  AQ  t  R(H)  then  by  theorem  4,  lim  | | 0  _  | f  *  since 
°  n^°  n~i 

lim  {B}n  I  *  0,  Observe  that  this  holds  irrespective  of  the  size  of 

n-H»  ° 


I AQl .  Now  the  error  in  the  iterates  is  obtained  as 


e  -  I  -  I  =  [B]  I  -  6  . 

n  o  n  o  n~l 


Note  that  the  norm  of  the  first  term  is  monotonically  decreasing  and 
thus  it  is  evident  that  the  algorithm  should  be  terminated  after  a 
certain  optimum  number  of  steps*  Unfortunately  the  exact  number  of 
iterations  depends  on  the  particular  Q  under  consideration  and  the 
growth  rate  of  I  l®n-iM  versus  the  decay  rate  of  ||[B]n 


